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1 Introduction 



Progress on our understanding of the neutrinos continues to be exhilarating. This progress is 
due mainly to experiments on neutrino oscillation. Here, we explain the physics of oscillation 
in vacuum and in matter. 



2 The physics of neutrino oscillation 



qq ■ Treatments of the physics of neutrino oscillation may be found in, for example, [1, 2]. Here, 
we give a slightly modified treatment, and explain some points that have caused puzzlement, 
such as the fact that, even if neutrinos are their own antiparticles, their interaction with 
matter can still cause a difference between neutrino and so-called "antineutrino" oscillations. 

We assume that the couplings of the neutrinos and charged leptons to the W boson are 
correctly described by the Standard Model, extended to take leptonic mixing into account. 
These couplings are then summarized by the Lagrangian 

c w = —% E (h a i x u^ Ll w; + vui x U* m t La W+) . (1) 

1=1,2,3 

Here, L denotes left-handed chiral projection, £ a is the charged-lepton mass eigenstate of 
flavor a (£ e is the electron, the muon, and £ T the tau), and Vi is a neutrino mass eigenstate. 
The constant g is the semiweak coupling constant, and U is the leptonic mixing matrix [3]. 
Supposing, as assumed by Eq. flTJ, that there are only three charged-lepton mass eigenstates, 
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and three neutrino mass eigenstates, U is 3 x 3, and may be written as 
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In the extended Standard Model, the 3x3 mixing matrix U is unitary, and we shall 
assume that this is also true in nature. However, we note that if there are "sterile" neutrinos 
(neutrinos that do not couple to the W or Z boson), then there are N > 3 neutrino mass 
eigenstates, and the leptonic mixing matrix U that is unitary is N x N, rather than 3x3. 
The 3x3 matrix of Eq. (|2J) is then just a submatrix, and is not unitary [4]. 

Supposing that the unitary mixing matrix is N x N, not because of the existence of sterile 
neutrinos but because there are N conventional lepton generations, how many physically- 
significant parameters does U contain? To see how many, we note first that an NxN complex 
matrix contains N 2 entries, each of which may have a real and an imaginary part. Thus, the 
matrix can be fully specified by 2N 2 real parameters. If the matrix is unitary, then each of its 
columns must be a vector of unit length: J2 a \U a i\ 2 — 1; i = 1,N. Together, these conditions 
are N constraints. In addition, each pair of columns in U must be orthogonal vectors: 
J2a U a j = 0; i,j — 1, N with i ^ j. Taking into account that each of these N(N — l)/2 
orthogonality conditions has both a real and an imaginary part, we see that these conditions 
impose N(N — 1) constraints. Thus, the number of independent parameters in a general 
N x N unitary matrix is 2N 2 — N — N(N — 1) = iV 2 . However, in the case of our unitary 
matrix, U, some of these parameters may be removed. From Eq. (JTj), (£ a \£wWiW~) oc U a i. 
Now, without affecting the physics, we are always free to redefine the state (£ a \ by multiplying 
it by a phase factor: (£ a \ — > (£' a \ = {£ c \e~ %LPa . Clearly, this has the effect of multiplying the 
U a i, f° r an h by the same factor: U ai — > U' ai = e~ tlfo, U a i. If there are N £ a , this phase 
redefinition of them may be used to remove N phases from U. It might be thought that 
analogous phase redefinition of the neutrinos z/j could be used to remove additional phases. 
However, unlike the quarks and charged leptons, the neutrino mass eigenstates z/j may be 
their own antiparticles: VI = z/j. This possibility motivates the search for neutrinoless nuclear 
double beta decay, as discussed at this school by K. Zuber. If VI = z^, then physically 
significant phases cannot be eliminated by phase redefinition of the z/j [5]. To allow for the 
possibility that V^ = z/j, we shall retain the phases that can be eliminated only when VI ^ z/j. 
Then U is left with iV 2 — iV physically significant parameters. These are commonly chosen 
to be "mixing angles" — parameters that would be present even if U were real — and complex 
phase factors. To see how many of the parameters are mixing angles, and how many are 
phases, let us imagine for a moment that U is real. Then it can be fully specified by its iV 2 
real entries. These are subject to the unitarity requirement that the N columns of U all have 
unit length: J2a U^i = 1, i = 1,N, and the requirement that all N(N — 1)/2 pairs of columns 
be orthogonal: J2 a U a iU a j = 0, i,j — 1, N with i ^ j. Hence, a real mixing matrix U for 
N generations has iV 2 — N — N(N — l)/2 = iV(iV — l)/2 physically significant parameters, 
and a complex one has this number of mixing angles. Since a complex U has iV(iV — 1) 
physically significant parameters in all, the fact that N(N — l)/2 of them are mixing angles 
means that the remaining N(N — l)/2 must be complex phase factors. 

In summary, a complex N x N unitary mixing matrix U for N lepton generations may 
contain — 



Neutrino Oscillation Phenomenology 



3 



N(N-l)/2 
N(N-l)/2 



mixing angles 
complex phase factors 



N(N-l) physically significant parameters in all 

Throughout most of these lecture notes, we will assume that N = 3. Then the mixing 
matrix contains three mixing angles and three complex phase factors. It can be shown that 
this matrix can be written in the form 



U 



1 

C23 S23 
— S23 C23 
e «i/2 
e^ 2 ' 2 
1 



x 



ci 3 s 13 e 
1 
a 



-si 3 e 



C13 



X 



C12 sia 
-S12 C12 
1 



(3) 



Here, cy = cos 6^- and Sjj = sinfly, where the are the three mixing angles. The quantities 
5, £1, and £ 2 are the three complex phases. 

From Eq. (pQ), we observe that the amplitude for the decay W + — > £ Q + z/j to yield the 
particular charged-lepton mass eigenstate £ Q in combination with the particular neutrino 
mass eigenstate is proportional to U*^ Thus, if we define the "neutrino state of flavor a", 
|i/ a ), with a = e, /x, or r, to be the neutrino state that accompanies the particular charged 
lepton £ a in leptonic W + decay, then we must have 



(4) 



1=1 



From Eq. ([T]), the amplitude for this v a to interact and produce the particular charged-lepton 
£p is proportional to 
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y9a 



(5) 



i=i 



where we have invoked the unitarity of U. We see that when a z/ e , the neutrino born in a 
W + decay that produced an e, interacts and produces a second charged lepton, the latter 
can only be an e. Similarly for and v T . 

We may invert Eq. (J3J) to obtain 



i) = E U aiWa) ■ 



(6) 



a=e,/i,T 



This expresses the mass eigenstate |z/j) in terms of the states of definite flavor, \u a ). We see 
that the flavor-a fraction of \vi) is simply |t/ Q j| 2 . 



2.1 Neutrino oscillation in vacuum 

Consider the vacuum neutrino oscillation experiment depicted schematically in the upper 
part of Figure dJ A neutrino source produces, via W exchange, the charged lepton £ a of 
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Figure 1. Neutrino flavor change (oscillation) in vacuum. "Amp" denotes an amplitude. 

flavor a, plus an accompanying neutrino that, by definition, must be a u a . The neutrino 
then propagates, in vacuum, a distance L to a target/detector. There, it interacts via W 
exchange and produces a second charged lepton lp of flavor (3. Thus, at the moment of its 
interaction in the detector, the neutrino is a Up. If the flavors a and f3 are different, then, 
during the neutrino's trip to the detector, it has changed, or "oscillated", from a u a into a 
v p . 

In the neutrino mass eigenstate basis, the particle that travels from the neutrino source 
to the detector is one or another of the mass eigenstates z/j. In a given event, we will not 
know which z/j was actually involved. Hence, the amplitude for the oscillation u a —>■ up, 
Amp [y a — > up), is a coherent sum over the contributions of all the Ui, as shown in the lower 
part of Figure HJ The contribution of an individual Ui is a product of three factors. The 
first is the amplitude for the neutrino produced together with the charged lepton £ a to be, 
in particular, a z/j. From Eq. this amplitude is U*^ as indicated in Figured) The second 
factor, Prop (z/j), is the amplitude for the mass eigenstate u^ to propagate from the source 
to the detector. The final factor is the amplitude for the charged lepton created when the 
Ui interacts in the detector to be, in particular, an lp. From Eq. (JT]), this amplitude is Upi. 

From elementary quantum mechanics, the propagation amplitude Prop(z/j) is simply 
exp [— irriiTi], where irti is the mass of Ui, and Tj is the proper time that elapses in the z/j 
rest frame during its propagation. By Lorentz invariance, m^Tj = Eit — piL, where L is the 
lab-frame distance between the neutrino source and the detector, t is the lab-frame time 
taken for the beam to traverse this distance, and Ei and pi are, respectively, the lab-frame 
energy and momentum of the z/j component of the beam. 

Once the absolute square |Amp {u a — > up)\ 2 is taken to compute the probability for the 
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oscillation v a — > up, only the relative phases of the propagation amplitudes Prop (z/j) for 
different mass eigenstates will have physical consequences. From the discussion above, the 
relative phase of Prop (vi) and Prop (i^), £0(12) is given by 

50(12) = (E 2 t-p 2 L)-(E 1 t- Pl L) 

= (p 1 -p 2 )L-(E 1 -E 2 )t . (7) 

In practice, experiments do not measure the transit time t. However, Lipkin has shown [6] 
that, to an excellent approximation, t may be taken to be L/v, where 

P1+P2 /o\ 

V = 



E 1 + E 2 

is an approximation to the average of the velocities of the v\ and v 2 components of the beam. 
We then have 

2 2 F 2 _ F 2 
5 0(12) ^-^L- 1 

Pi + P2 Pi + P2 

= (m\-m 2 ) — - — = {ml-m\)-^- , (9) 
P1+P2 2E ' v ' 

where, in the last step, we have used the fact that for highly relativistic neutrinos, pi and 
p 2 are both approximately equal to the beam energy E. We conclude that all the relative 
phases in Amp( v a — > vp) will be correct if we take 

Prop(z/ 4 ) = e - im " L/2E . (10) 

Combining the factors that appear in the lower part of Figure [U we have 

Amy J (v a ^v p ) = Y,U: i e- m * L / 2E Up l . (11) 

i 

Squaring, and making judicious use of the unitarity of U, we find that the probability of 
v a -> up, P (u a -> up), is given by 

P(u a -> up) = I Amp (u a -> up)] 2 

= 5 a p - 4 ]T Re (U^U^U^) sin 2 ! (Am 2 .L / 'AE) 

i>j 

+ 2j2l^(U* m U Pl U a3 U; 3 )sm(Am 2 J L/2E) . (12) 

i>j 

Here, Am^ = mf — m 2 is the splitting between the squared masses of and Uj. It is clear 
from the derivation of Eq. ffT2"]) that this expression would hold for any number of flavors 
and equal number of mass eigenstates. 

Given that the particles described by the oscillation probability of Eq. (fT2l are born with 
an £ a and convert into an lp in the detector, they are neutrinos, rather than antineutrinos 
(should there be a difference). To obtain the corresponding oscillation probability for an- 
tineutrinos, we observe that TT^ — > T7p~ is the CPT-mirror image of vp — > v a . Thus, if CPT 
invariance holds, 

P(^ - vf) = P(vp -> v a ) . (13) 
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Now, from Eq. (fT2|) . we see that 

P(f/j -> ^; C/) = PK -> ^; C/*) • (14) 

Hence, assuming CPT invariance holds, 

Pfc -+T7p-U) = P(v a -> ^; £/*) . (15) 

That is, the probability for oscillation of an antineutrino is the same as that for a neutrino, 
except that the mixing matrix U is replaced by its complex conjugate. Thus, from Eq. (TT2]) . 

P(fc2-^) = ^-4^Re(^^[/ aj [/* i )sin 2 (Am2.L/4 J E;) 

( - , 2 ^ m (UaiUpiUajUpj) sin(Am 2 L/2£) . (16) 

We see that if U is not real, the probabilities for v a — > i/^ and for the corresponding an- 
tineutrino oscillation, z/^ — > I7g, will in general differ. Since z/ a — > vp and — > X/g are 
CP-mirror-image processes, this difference will be a violation of CP invariance. 

As Eq. (1161) makes clear, neutrino oscillation in vacuum from one flavor a into a different 
one j3 implies nonzero mass splittings Am|, hence nonzero neutrino masses. It also implies 
nontrivial leptonic mixing. That is, the mixing matrix U cannot be diagonal. 



Including the so-far omitted factors of ft and c, we have 

A K ± = L^CV)^ . (17) 

From Eq. ffTBT) . if the U matrix cooperates, the probability for u a — > up, ^ a, will be 
appreciable if the kinematical phase difference in Eq. (fTTl) is 0(1) or larger. This requires 
only that for some ij, 

Thus, for example, an experiment that studies 1 GeV neutrinos that travel a distance 
L ~ 10 4 km, the diameter of the earth, will be sensitive to neutrino (mass) 2 splittings Amfj 
as small as 10~ 4 eV 2 . Through quantum interference between neutrino mass eigenstates of 
different masses, neutrino oscillation gives us sensitivity to very tiny (mass) 2 splittings. How- 
ever, as Eq. (fIBl) underscores, oscillation cannot determine the masses of the individual 
mass eigenstates. To learn those will require another approach. 

There are basically two kinds of neutrino oscillation experiments. In the first, an ap- 
pearance experiment, one starts with a beam of neutrinos that initially are purely of flavor 
a, and looks for the appearance in this beam of neutrinos of a new flavor (3, (3 ^ a, that 
were not originally present in the beam. In the second kind of experiment, a disappearance 
experiment, one starts with a known flux of u a , and looks to see whether some of the initial 
u a flux disappears as the beam travels. 

By the definition of "probability", the probability that a neutrino changes flavor, plus 
the probability that it does not change flavor, must equal unity. That is, we must have 

^P(z/ a ^z/ /3 ) = ^P(x^^I^) = l , (19) 
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where the sum is over all final flavors j3, including the initial flavor a. From the unitarity of U, 
which implies that Y,p UpiUp, = 5^-, it immediately follows that the oscillation probabilities 
of Eq. ([TBI) do obey this constraint. 

Neutrino flavor oscillation does not change the total flux in a neutrino beam. It merely 
redistributes it among the flavors. However, if we create a beam of neutrinos that at birth 
are of some active (i.e., weakly interacting) flavor, (muon neutrinos, for example), and some 
of these neutrinos oscillate into sterile (i.e., non- interacting) flavors, then some of the total 
active neutrino flux will have disappeared. 

The combination of the CPT-invariance constraint of Eq. (fl3|) and the probability con- 
straint of Eq. (Tl9l has powerful consequences for CP violation. To see this, consider the 
CP-violating differences 

A a p = P(u a -> up) - P(T^ -> Up) . (20) 
If CPT invariance holds, then from Eq. ffl3l) 

Ap a = -A a p . (21) 

In particular, 

A aa = . (22) 

That is, there can be no CP- violating difference between the survival probabilities P(u a — ► 
u a ) and P(Pq — > u^)- Hence, there can be no observable CP violation in a disappearance 
experiment. Now, from Eq. ( fl9l) . it follows that 

E A ^ = ° > ( 23 ) 

P 

where the sum runs over all flavors, including (3 = a. However, in view of Eq. (T22]) . Eq. ff23l) 
implies that 

E A ^ = ° • (24) 

P^a 

If there are only three neutrino flavors, u e , u^, and u T , then this constraint implies that, in 
particular, 

A eM + A er = and A Me + A MT = . (25) 

From these relations and Eq. (l2Tj) . we see that 

A efl = A Mr = A Te = -A Me = -A rM = -A er = A . (26) 

In summary, if CPT holds, then the CP-violating difference A a p = P(u a — * up) — P(v^ — > up) 
can be nonvanishing only for (3 ^ a. If, in addition, there are only three flavors, then the 
six possibly-nonvanishing A a p, shown in Eq. fT26|) . must all be equal, apart from a predicted 
minus sign [7]. (If there are more than three flavors, then Eq. fT26|) need not hold.) 

Counter to intuition, the CP- violating difference A a p = P(u a — > up) — P(Pq — > up) 
between neutrino and what we conventionally call "antineutrino" oscillation probabilities 
can still be nonvanishing even when the u\ are identical to their antiparticles. Indeed, A a p 
is actually completely independent of whether the Ui are their own antiparticles or not. We 
illustrate this by comparing the processes u^ — > u e and "T^ — > u^" , depicted in Figure [2j 
In Up — > u e , the neutrino is created together with a fi + in tt + decay. After traveling down 
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Figure 2. T/ie CP-mirror-image oscillations v u 



and 



V7, 



u e ". In each process, 



the particle that travels down the beamline is one or another of the mass eigenstates, and 
the amplitude is a coherent sum over the contributions of these eigenstates, as indicated. In 
"u^ — > V e " , the mass eigenstate VI may or may not be identical, apart from its polarization, 
to the corresponding v± in — > v e . The propagator for this VI, ex~p(—im?L/2E), is identical 
to that for the corresponding in either case. The elements of the U matrix that, according 
to Eq. (QJ) ; appear at the beam-particle production and detection vertices are shown. 



a beamline to a detector, it is detected via its production of an e~. In the corresponding 
"antineutrino" oscillation, "T^J — > Vf , the particle that travels down the beamline is created 
together with a pT in n~ decay, and is detected via its production in the detector of an e + . 
One never directly observes the particle that travels down the beamline; it is an intermediate 
state. In terms of the charged leptons that one does ( or at least can) observe, — > v e and 
U V^ — > V~z are clearly different, CP-mirror-image processes: the first involves a fi + and e~, 
while the second involves a fi~ and e + . Thus, even if V\ = z/j, — > v e and "PJJ — > V^" can 
have different probabilities, and if they do, the difference is a violation of CP invariance. 

Even if V\ — z^, the beam particle will create an e~ in —>■ u e , but an e + in U V^ — > v£\ 
because it is oppositely polarized in the two processes. Due to the chirally left-handed 
structure of the weak interaction, reflected in the Lagrangian of Eq. (TjQ), the beam particle 
will have helicity h = —1/2 in the first process, but h = +1/2 in the second. Due to this 
same parity- violating left-handed structure, the h = —1/2 beam particle will create an e~ 
(via the first term in Eq. (TjQ)) in — > v e , while the h = +1/2 beam particle will create an 
e + (via the second term in Eq. (CQ)) in U V^ —>■ V^\ 

From the amplitude factors displayed in Figure [U we see that while 

Km V {v^v e )=Y^U; i e- m ^U et , (27) 
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Amp(^ 



(2* 



These expressions hold whether = v$ or not. Thus, in either case, if the CP-violating phase 
8 in Eq. ([3]) is not zero or n, so that U is complex, the interference terms in P(z/ M — > u e ) and 



will differ. As a result, the CP-violating difference P(f (U 



will 

be nonzero. Furthermore, the value of this difference will not depend on whether V[ = Ui, and 
this value will be correctly implied by Eq. (fT6l) . which holds regardless of whether = Vi. 

The general expression for P(fe ) — > Vg), Eq. ( fl6l) . simplifies considerably in some impor- 
tant special cases. One such case is the simplified world in which there are only two charged 
leptons, say e and /x, two corresponding neutrinos of definite flavor, v e and z/ M , and two 
neutrino mass eigenstates, V\ and i>2, that make up v e and v^. From our earlier analysis of 
the number of parameters in a mixing matrix, we know that the 2x2 unitary mixing matrix 
U for this two-flavor world may contain one mixing angle and one complex phase factor. It 
may easily be shown that U may be written in the form 



u = 




u e2 ' 




cos 9 sin 9 




" e^ 2 


" 






X 










U, 2 _ 




— sin 9 cos 9 







I 



(29) 



where 9 is the mixing angle and £ is the phase. With Am^ = Am 2 the sole (mass) 2 splitting 
in the problem, we find from the U of Eq. (|29|) and the general expression of Eq. (fl6|) that 



P(W -> ^) = P(^ -> = sin 2 2#sin 2 (Am 2 L/4£) 



(30) 



and that 



P(fc£ -> fc?) = P(i^ -> ^) = 1 - sin 2 2#sin 2 (Am 2 L/4£) 



(31) 



As we know, the real world contains (at least) three charged leptons £ a , three corre- 
sponding neutrinos of definite flavor u a , and three underlying neutrino mass eigenstates Ui 
that make up the v a . Thus, the two-neutrino oscillation formulae of Eqs. ( 1301) and (l3~Pj) do 
not apply. However, if there are only three flavors, then under certain circumstances, rather 
similar simple formulae do apply. To see this, we note that the three-neutrino (mass) 2 spec- 
trum has been observed to have the form shown in Figure [3] [2]. The splitting Am^, which 
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\ V- 






i 




(mass) 2 




|Arr 


2 32 l or 




\= 




* ^Am 2 2l 



|Am 2 



311 



■Am 



21 



Figure 3. The three-neutrino (mass) 2 spectrum. 



Am 2 ^ which 



drives the behavior of solar neutrinos, is roughly 30 times smaller than Am 2 2 
drives the behavior of atmospheric neutrinos. (It is not known whether the closely-spaced 
pair V1-V2 is at the bottom or the top of the spectrum.) If an experiment is performed with 
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L/E such that Am^ 2 L/E = 0{l), then Am^ L/E <C 1, and in first approximation, this ex- 
periment cannot "see" the small splitting Am^. Neglecting this small splitting in Eq. (fTBj) . 
this equation and the unitarity of U imply that, for j3 ^ a, 

P(^^^)^4\U a3 U p3 \ 2 S m 2 (Am 2 2 L/4E) . (32) 

Similarly, they imply that, for /3 — a, 

P(fr2^tt) = l-MU a3 \ 2 (l-\U a3 \ 2 )sm 2 (Aml 2 L/4E) . (33) 

We see that, by measuring these simple oscillation probabilities, experiments with Am^ L/4E 
= 0(1) can determine the flavor content of the isolated member of the spectrum, u 3 . 

2.2 Neutrino oscillation in matter 

Inside matter, the coherent forward scattering of neutrinos from the electrons, protons, and 
neutrons that make up the matter leads to neutrino effective masses and mixing angles 
that differ from their vacuum counterparts. As a result, inside matter, the probabilities for 
neutrino oscillations differ from their vacuum counterparts. 

The Standard-Model interactions between neutrinos and other particles do not change 
flavor. Thus, barring hypothetical non-Standard-Model flavor- changing interactions, the 
observation of neutrino flavor change implies neutrino mass and leptonic mixing, even if the 
observation involves neutrinos passing through matter. 

Neutrino propagation in matter may be conveniently treated via the laboratory-frame 
Schrodinger time-evolution equation 

i^*(t)=W*(t) . (34) 

Here, t is the time, and is a multi-component neutrino wave function. Its a component, 
ty a (t), is the amplitude for the neutrino to have flavor a at time t. If there are N flavors, 
the Hamiltonian TC is an N x N matrix in flavor space. In matter, this matrix includes 
interaction energies arising from neutrino-matter interactions mediated by W or Z exchange. 
According to the Standard Model, the Z-mediated interactions neither change neutrino flavor 
nor depend on the flavor. Thus, they add to TC a term proportional to the identity matrix. 
Such a term shifts all the eigenvalues of TC by a common amount, leaving the splittings 
between the eigenvalues unchanged. Now, as we have seen when discussing neutrino flavor 
oscillation in vacuum, the amplitude for oscillation depends only on the relative phases of the 
different neutrino eigenstates. This means that it depends only on the splittings between the 
eigenvalues, and will not be affected by an interaction that merely shifts all the eigenvalues 
by the same amount. Thus, if our purpose is to treat neutrino flavor oscillation, we may 
omit the Z-exchange contribution to TC. 

The VF-exchange contribution is another matter. From the Standard Model, it follows 
that coherent forward i/ e -electron scattering via the M^-exchange diagram of Figure H] adds 
to the v e -v e element of TC, TC UeUe , an interaction energy 

V = V2G F N e . (35) 

Here, Gf is the Fermi coupling constant, and N e is the number of electrons per unit volume 
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Figure 4. The W -exchange interaction that modifies neutrino flavor oscillation in matter. 



in the matter through which the neutrinos are passing. The Fermi constant appears in 
V because it is a measure of the amplitude for the diagram in Figure HI and the density 
N e appears because the coherent scattering amplitude will obviously depend on how many 
electrons are present to contribute. The Standard Model tells us that for antineutrinos in 
matter, V is replaced by —V. 

Since v e is the only neutrino flavor that couples to an electron and a W, IF-mediated 
v — e scattering affects only the v e -v e element of 7i. Thus, its contribution to 7i is not 
proportional to the identity matrix, and it does affect neutrino flavor oscillation. 

Neutrino flavor change in matter is illustrated by the case where there are only two 
significant flavors, say v e and and, correspondingly, two significant mass eigenstates. 
The Hamiltonian 7i in the Schrodinger equation, Eq. is then a 2 x 2 matrix in v e -v^ 

space. Taking into account v e -e scattering via W exchange, but omitting some irrelevant 
contributions that are proportional to the identity matrix, we readily find [1] that 



H 



AE 



cos 29 m 
sm29 M 



sin 29 m 
cos 20 M 



(36) 



Here, E is the energy of the neutrinos, and Am 2 M and 9m are, respectively, the effective 
(mass) 2 splitting and the effective mixing angle in matter. These effective quantities are 
related to their vacuum counterparts, Am 2 and 9, by 



and 



In these expressions, 



Am 2 M 



sin 2 29 



Am V sin 2 29 + (cos 29 - x v 



M 



sin 2 29 
sin 2 29 + (cos 29 



X j, 



2V2G F N e E 
Am 2 



(37) 
(38) 

(39) 



is a measure of the importance of the effects of matter. In vacuum, N e and consequently 
x u vanishes, and, as confirmed by Eqs. (1371) and (1381) . Am 2 M and 9m revert to the vacuum 
values, Am 2 and 9, respectively. 

Imagine an accelerator-generated neutrino beam that travels a distance L ~1000km 
through the earth to a detector. The electron density iV e encountered by this beam will 
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be that of the earth's mantle, and approximately constant. Then x v , Am 2 M , 6m, and 7i 
will all be ~ position- independent. From Eqs. and and straightforward quantum 
mechanics, it follows that 

P(^ e -> ^) = PK ^) = sin2 sin 2 (Am^L/4E7) . (40) 

This is the usual two-neutrino oscillation result, Eq. (I5U|) . except that the vacuum parameters 
6 and Am 2 are replaced by their counterparts in matter, 6 m and Am M . If N e — > 0, so that 
x u — > 0, the oscillation probabilities in matter of Eq. ([40]) become the vacuum probabilities 
of Eq. (1301) . as they must. 

The size of the effect of matter may be judged by the size of x v . For the illustrative 
beam that we are considering, the actual three-neutrino vacuum (mass) 2 splitting Am 2 that 
will most strongly influence flavor oscillation is probably the large one, Am^ = Am 2 2 . 
Experimentally, Am^ ~ 2.4 x 10~ 3 eV 2 [8]. For this Am 2 , we find from Eq. fl39l) that 

\x u \ ~ £/12GeV . (41) 

Thus, for E = 0.5 GeV, the matter effect is quite small, for E = 2 GeV it is modest, and for 
E = 20 GeV it is large. 

As already mentioned, when antineutrinos, rather than neutrinos, propagate through 
matter, the interaction energy V is replaced by — V. It follows readily that, as a result, x v , 
Eq. ([39]), is replaced in Eqs. (|57j)-(P]l by 

Xjj = Xy . (^^) 

We see that this change has the consequence that, within matter, the effective (mass) 2 
splitting and the effective mixing angle for antineutrinos are different than they are for 
neutrinos. As a result, within matter the flavor oscillation of an antineutrino beam will differ 
from that of a neutrino beam. The two-flavor Hamiltonian TC is still given by Eq. (155]) . and 
the two-flavor oscillation probability in matter of constant density is still given by Eq. f)40p . 
but the quantities Am 2 M and 6m have different values than they did in the neutrino case. 

Earlier, we raised the possibility that neutrinos are their own antiparticles. That is, we 
imagined that, for a given momentum p and helicity h, perhaps each neutrino mass eigenstate 
Vi is identical to its antiparticle: Ui(p,h) = Vj(p,h). Suppose that this is indeed the case. 
Will the interaction with matter still cause the flavor oscillation of an "antineutrino" beam 
to differ from that of a "neutrino" beam within matter? In practical terms the answer is 
"yes". The reason is that, in practice, the "neutrino" and "antineutrino" beams that we 
study are never of the same helicity. A "neutrino" is the particle V" produced, for example, 
in the decay W + — > e + + v. As already noted, owing to the chirally left-handed structure 
of the weak interaction, this 'V will be of left-handed helicity: h = —1/2. In contrast, an 
"antineutrino" is the particle "z/" produced in W~ — > e~ + v. As already noted, owing to the 
same structure of the weak interaction, this u v n will be of right-handed helicity: h = +1/2. 
Since the weak interaction is not invariant under parity, the interaction in matter of the 
left-handed 'V and the right-handed "z/" will be quite different, even if helicity is the only 
difference between the V" and the "z/" . Only the first term on the right-hand side of Eq. (TjQ) 
can couple an incoming left-handed beam particle to an electron, while only the second 
term can couple an incoming right-handed beam particle. These two terms lead to different 
scattering amplitudes. These amplitudes do not depend on whether the "z/" and "z>" beams 
differ only in helicity, or in some other way as well. 
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Future accelerator neutrino experiments hope to study — > z/ e and — > I^T in matter 
under conditions where all three of the known neutrino mass eigenstates ^1,2,3, or equiva- 
lently both of the known splittings Am^ and Am| 1; play significant roles. The oscillation 
probabilities are then more complicated than the expression of Eq. (j4"01) . However, since 
a = Am^/Am^ ~ 1/30 [2] and sin 2 2^13 < 0.2 [9], the probability for z/^ — > v e in matter is 
well approximated by [10] 

PiyvL -> ^e) = sin 2 2fl 13 Ti - a sin 2# 13 T 2 + a sin 2fl 13 T 3 + a 2 T 4 , (43) 

where 

_ .2 fl sin 2 [(l-q;,)A] 

" sm (1 _^ )2 , (44) 

T 2 = sin J sin 26 12 sin 26 23 sin A sin ^ A ) sin p Z ^ , (45) 

Xy (1 ^//) 

m sin(xjyA) sin[(l — x„)Al . . 

T 3 = cos 5 sin 2^ 12 sin 2^3 cos A — K -0—L — l f- ^—^ , (46) 



and 



T 4 = cos 2 g 2 3sin 2 2g 12 Sin2( ^ A) . (47) 

XT, 



In these expressions, A = Am^L/AE is the kinematical phase of the oscillation, and x v is 
the matter-effect quantity defined by Eq. (1391) . with Am 2 now taken to be Amf^ In the 
appearance probability P(^ — > f e ), the Ti term represents the oscillation due to the splitting 
Am 2 1; the T4 term represents the oscillation due to the splitting Am 21 , and the T 2 and T 3 
terms are the CP-violating and CP-conserving interference terms, respectively. 

The probability for the corresponding antineutrino oscillation, P(T^ — > T^), is the same 
as the probability P(^ — > v e ) given by Eqs. P3l - fT4Tj) . but with x u replaced by x„ = —x v 
and sin S by — sin 5: both the matter effect and CP violation lead to a difference between 
the — > v e and —>■ V~ e oscillation probabilities. In view of the dependence of x v on Am^, 
and in particular on the sign of Am 2 1; the matter effect can reveal whether the neutrino 
mass spectrum has the closely-spaced v x -v 2 pair at the bottom or the top (see Figure [3]). 
However, to determine the nature of the spectrum, and to establish the presence of CP 
violation, it obviously will be necessary to disentangle the matter effect from CP violation in 
the neutrino-antineutrino oscillation probability difference that is actually observed. To this 
end, complementary measurements will be extremely important. These can take advantage 
of the differing dependences on the matter effect and on CP violation in P(f M — > v e ) and 
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